Optical solitons in PT-symmetric nonlinear couplers with gain and loss 
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We study spatial and temporal solitons in the VT symmetric coupler with gain in one waveguide 
and loss in the other. Stability properties of the high- and low-frequency solitons are found to 
be completely determined by a single combination of the soliton's amplitude and the gain/loss 
coefficient of the waveguides. The unstable perturbations of the high-frequency soliton break the 
symmetry between its active and lossy components which results in a blowup of the soliton or a 
formation of a long-lived breather state. The unstable perturbations of the low-frequency soliton 
separate its two components in space blocking the power drainage of the active component and 
cutting the power supply to the lossy one. Eventually this also leads to the blowup or breathing. 



PACS numbers: 42.25.Bs, 11.30.Er, 42.82.Et 

I. INTRODUCTION 

Optical solitons are formed when nonlinear effects com- 
pensate the diffractive broadening of light beams (spatial 
solitons) or dispersive spreading of optical pulses (tem- 
poral solitons). Although these localized structures arise 
both in conservative settings and in systems with active 
and lossy elements, properties of dissipative optical soli- 
tons [l] [2] show significant differences from those of their 
conservative counterparts [3 . In particular, the ampli- 
tudes of solitons in conservative systems are free to vary 
over continuous ranges whereas the generic dissipative 
systems can only support solitons at special amplitude 
values determined by the balance between gain and loss. 

In an interesting turn of events, it was realised recently 
that there is a class of optical systems where dissipative 
solitons arise in continuous families. These systems con- 
sist of elements with gain and loss arranged in a par- 
ticular symmetric way [4 . The symmetry here can be 
interpreted as an optics equivalent O |6] of the VT 
(parity-time) symmetry in quantum mechanics [THTQ]. 
The T^T-symmetric potentials in quantum mechanics are 
essentially complex potentials which however exhibit a 
purely real spectrum of energies, with the implication 
that their time-dependent eigenfunctions show no decay 
or growth. Despite this similarity to the hermitian quan- 
tum mechanics, the T^T- symmetric quantum systems dis- 
play a variety of anomalous phenomena stemming from 
the nonhermitian mode interference [11-16 . 

The first experimental demonstrations of the VT- 
symmetric effects in optics were in two-waveguide di- 
rectional linear couplers composed of waveguides with 
gain and loss [17^18^. Theoretical analyses suggest that 
such couplers, operating in the nonlinear regime, can be 
used for the all-optical signal control p!9H22] . Arrays of 
the PT-symmetric couplers were proposed as a feasible 
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FIG. 1. A schematic of PT-symmetric coupled waveguides 
with gain (top waveguide) and loss (bottom waveguide), 
(a) Two waveguides on the plane where t denotes the logitu- 
dinal and x the transversal spatial coordinate. The transverse 
x-profiles of the stationary light beams evolve as they extend 
in the t direction, (b) A pair of one-dimensional waveguides 
where light pulses undergo temporal evolution as they travel 
along the x axis. 



means of control of the spatial beam dynamics, including 
the formation and switching of spatial solitons [23H25] . 

This paper is concerned with the T^T-symmetric cou- 
plers with an extra spatial or temporal degree of freedom. 
We consider the situation of stationary light beams in the 
coupled planar waveguides [i.e. waveguides extended in 
the transverse direction, Fig.jlja)] and that of the optical 
pulses in coupled one-dimensional waveg uides [Fig.jljb)]. 
The configuration shown in Fig. [TJa) can also be seen as 
the strong-coupling limit of the array of coupled dimers 
discussed in the recent Ref. [25 . 

Our study focuses on spatial and temporal solitons in 
such couplers. The reader should be alerted up front 
that we use the term "soliton" simply as a synonym for 
solitary wave or localised pulse. No a priori stability is 
implied by the use of this term. It is the objective of this 
study to classify the T^T-symmetric solitons into stable 
and unstable ones. 

In addition to the analysis of the solitons' stability 
properties and numerical study of the linearisation eigen- 
values, we uncover the instability mechanisms and simu- 
late the nonlinear evolution of the unstable solitons. 

The paper is organized as follows. In the next sec- 
tion (Sec. |ll|) we formulate the mathematical model and 
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identify physically meaningful integral characteristics of 
the associated evolution. Two families of high- and low- 
frequency bright soliton solutions of the model are intro- 
duced in Sec. 



Whether Eqs.(2.1) are employed to describe the sta- 
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In the following section, Sec. |IV[ we 
outline the general framework for their stability analysis. 
The stability eigenvalues of the high-frequency soliton are 
classified in section |V| there, we also follow the nonlin- 
ear evolution of instability when the soliton is found to 
be unstable. Section |VI] contains a similar study of the 
low-frequency soliton. 



Finally, section |VII| summarises 
conclusions of our work while three Appendices detail 
mathematical analyses of the stability eigenvalues. 



II. THE MODEL 

To describe the dynamics of stationary light beams 
and pulses in coupled waveguides illustrated in Fig. [l] 
we extend the equations of the nonlinear T^T-symmetric 
coupler [I9H2T]. In the physical setting of Fig. [ija), our 
extension takes into account the diffraction of the beam 
while in the situation represented by Fig.jljb), we extend 
the system to include the effect of the pulse dispersion. 
The resulting equations have the following dimensionless 
form: 

iut + Uxx + 2\u\'^u = —V -\- iju^ 
ivt + Vxx + 2\vfv = -u- i'jv. (2.1) 

We note that this model was originally introduced as 
the continuous limit for a one-dimensional array of VT- 
symmetric cou pled waveguides [25] , 

In Eqs.( |2.l| , the u and v variables are the normal- 
ized complex mode amplitudes in the top and bottom 
waveguides of Fig. [l] Considering stationary beams in 
the planar geometry of Fig. [ija), the t variable is the 
(spatial) coordinate in the propagation direction while x 
is the transversal coordinate. In the temporal pulse in- 
terpretation [Fig. [ijb)], t stands for time and x for the 
spatial coordinate in the frame moving with the pulse 
group velocity. Here we assume that the group veloci- 
ties and second-order dispersions in the waveguides are 
matched so as to satisfy the T^T-symmetry condition. 
By scaling the x variable properly, we have normalized 
the coefficients in front of Uxx and Vxx terms to unity. 
(These terms account for the diffraction of spatial beams 
and dispersion of temporal pulses.) 

We assume that the Kerr nonlinearity coefficients in 
the two waveguides (coefficients in front of \u\'^u and 
\v\'^v) have the same value as this is necessary for the 
existence of solitons [25 . (These coefficients have been 
normalized to 2 by scaling the mode amplitudes.) To en- 
sure the existence of bright solitons, we take equal signs 
in front of the diffraction/dispersion and nonlinear terms. 
This corresponds to the self-focusing nonlinearity in the 
case of beams [Fig.[TJa)], and to the anomalous dispersion 
with positive Kerr nonlinearity or normal dispersion with 
negative nonlinearity — in the case of pulses [Fig. [TJb)]. 



tionary planar beams or temporal pulses, the first terms 
in the right-hand sides of (2.1) account for the coupling 



between the modes propagating in the two waveguides. 
The 7-terms describe the gain in one and loss in the other 
waveguide. Without loss of generality 7 can be taken pos- 
itive; this choice corresponds to the gain in the top and 
loss in the bottom waveguide. The gain and loss coeffi- 
cients are taken equal to conform to the T^T-symmetry 
condition [18 . 

We close this section by noting several physically 
meaningful quantities which prove useful in the under- 
standing of the dynamics described by Eqs.(2.1). The 



first pair of variables give the powers associated with the 
u and V fields, respectively: 



I 



dx. 



(2.2) 



Neither individual powers nor their sum are conserved if 
7 7^ 0; however, the rate of change of the total power has 
a simple and insightful expression: 

^(K+n) =27(K-n). (2.3) 
The momenta carried by the u and v components, 

Mu = ^ J {ulu-Uj:U')dx, Mv = ^ J {Vx^-Vxv'^)dx, 

(2.4) 

are not conserved either. The total momentum satisfies 
d 



dt 



{Mu + Mu) = 2j{Mu-Mv). (2.5) 



Finally we note the rate equation 
dH 



dt 



27(7^„-7^„), 



(2.6) 



where 



'H = I + - {\u\* + \v\^) - {vu* + v*u)] dx 



and 



Hu = 



2\u\'^)dx, = 



v = J{\Vx 



(2.7) 



2\v\'^)dx. 



(2.8) 

The integral H plays the role of the Hamiltonian of 
Eqs.(2.1) in the situation where there is no loss or gain 

(7 =or 



An immediate consequence ofEqs.(2.3), (2.5) and (2.6) 



is that all stationary states in the system (2.1) have to 
display symmetry between their two components: Vu = 



III. SOLITONS 

Pr oceeding to the analysis of solutions to the system 
(2.1), it is convenient to make a change of variables 



u{x,t) = e'^^'-^^U{x,t), v{x,t) = e'^'V{x,t), (3.1) 
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where ^ is a constant angle satisfying 
sin 6 — 7, 



(3.2) 



and O is an arbitrary real parameter which will be con- 
veniently chosen later. The transformation (3.1) casts 
equations (2.1) in the form 



iUt + U:cx -nU + 2\U\^U = -cos0V + ij{U - V), 

iVt + V:cx-^V + 2\VfV = -cos eU + ij{U-V). 

(3.3) 



The system (3.3) admits an obvious reduction U = 
V = <p to the scalar cubic Schrodinger equation, 



i(pt + (pxx-a'^(p + ^\<P\^<P = 0, 



(3.4) 



where = Q — cosO. The relation (3.2) defines two 
different angles, 6 = arcsin7 and = n — arcsin7. Ac- 



cordingly, Eq.(3.4) d escri bes two separate invariant man- 
ifolds of the system ( |2.1[ ). Both invariant manifolds are 
characterise d by the Hamiltonian evolution. 

Equation (3.4) has a family of stationary soliton solu- 
tions. Without loss of generality we can restrict ourselves 
to time-independent solutions, 

<p{x) = a sech(ax). 

These define two coexisting families of stationary solitons 
of the original system (2.1), with arbitrary amplitudes 



a > and the corresponding frequencies 
ft = -\- cos 0. 

The two families of solitons are distinguished by the 
sign of cos^. One has cos^ = y^l — 7^ > 0; we will be 
denoting the corresponding solitons by = (11+, 
The other one [denoted i/j- = {u-,V-) in what follows] 
has cos^ = — Y^l — 7^ < 0. Note that for the given 
amplitude a, the frequency Q corresponding to the soliton 
is higher than that for ip^ . For this reason, we will be 
referring to the two solitons as 'high frequency' and 'low 
frequency' solitons. (It is fitting to note that the previous 
authors [25] considered the high-frequency soliton only.) 

We note that either family exists only if 



7 < 1. 



(3.5) 



At the same time, Eq.(3.5) gives the stability condition 
for the background solution u = v = 0. In what follows, 
we will assume that (3.5) is always imposed. 



IV. STABILITY FRAMEWORK 

A. Perturbation decomposition 

To classify the stability of the two families of solitons, 
we let 

U{x, t) = 0(x) + SU{x, t), V{x, t) = 0(x) + SV{x, t) 

(4.1) 



and linearise Eqs.(3.3) in 5U and 5V. As we will see, a 



special role is played by the symmetric and antisymmet- 
ric combinations 



P = 



5 U + 5V 



5U-5V 



(4.2) 



Since the linearised equations are autonomous in time, 
it is sufficient to consider separable solutions of the form 



P ^ 



^[(Pi + ^^2) coscjt + {pi + ip2) sinut] 



Q = ^""^[{Qi + ^<72) coscjt + ((7i + iq^) sincjt], (4.3) 

where we have introduced real components of four com- 
plex functions pi{x), P2{x), qi{x) and q2{x): 



Pi -- 
Qi 



-Pi - 



P2 
(12 



P2 



■W2^ 



In (4.3), both u and uj are assumed to be real. 



Substituting (4.3) in the linearised equations yields an 



eigenvalue problem 

{C — cos 6)p + 277^ = /i Jp, 
{C + cos 0)q = jiJq 

for two-component complex vectors 



(4.4a) 
(4.4b) 



P 



In ( |4.4a[ )-( |4.4b| ) we have defined ji = v — iuj and intro- 



duced the operator 
C = 



d^/dx^^n- 






-d^/dx'^ ^n-2(i)^ 



The J in the right-hand sides of ( |4.4a )-(4.4b) stands for 
a skew-symmetric matrix: 



The rotation (4.2) did not diagonalise the 2x2 block 



supermatrix in the left-hand side of (4.4); however it 



brought it to the triangular form. The triangular block 
matrix 



-cos^ 27 J 
£ + cos ^ 



has two eigenvectors. One is 



with p satisfying 



{jC — cos 6)p = 11 Jp. 



(4.5) 



(4.6) 



(4.7) 



This is nothing but the linearised eigenvalue problem for 
the unperturbed cubic Schrodinger equation (the inte- 
grable nonlinear Schrodinger equation). (Thus p is the 
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component of the perturbation that hes in the tangent 
space to the conservative manifold containing the soh- 
ton.) The spectrum of ja consists of a four- fold zero 
eigenvalue and the continuous spectrum which lies on 
the imaginary axis; there are no unstable /i's here. 
The second eigenvector, 



(2.8), the right-hand sides of (2.3), (2.5) and (2.6) are 



(4.8) 



includes a nonzero q component. These q arise as eigen- 
functions of the operator (|4.4b[). When ^ ^ 0, Eq.(|4.4a) 



becomes a nonhomogeneous equation with the right-hand 
side determined by q: 

mp= -2-fJq. (4.9) 

Here is a nonsymmetric operator defined by 

or = £ — cos — II J. 

As discussed in the previous paragraph, the eigenvalue 
problem (4.7) does not have nonzero eigenvalues; hence 



the adjoint operator 



C — cosO + /iJ 



has null ei genv ectors only if /i = 0. Therefore, for any 
/i 7^ 0, Eq.(4.9) has a bounded solution. 



Thus the stability analysis reduces to solving the eigen- 
value problem (4.4b). It is important to emphasise that. 



despi te the prese nce o f gain and loss, the eigenvalue prob- 
lems (4.4b) and (4.7) are symplectic (that is, pertaining 
to Hamiltonian evolutions). The implication is that sta- 
ble limit cycles are not the objects that can be expected 
to bifurcate from stationary solitons when the latter lose 
their stability. That is, the instability should evolve ac- 
cording to scenarios that are characteristic for conserva- 
tive systems, e.g. breakup into long-lived transient struc- 
tures, singular growth etc. 

It is also worth commenting on the significance of the 
triangular representation (4.4) which decomposes small 



perturbations into the part tangent to the conservative 
invariant manifold containing the soliton, and the part 
that is transversal to this manifold. This decomposition 
alone explains some numerical observations of the previ- 
ous authors [25] , in particular the instability of the bond- 
centred solitons in the discrete case. Indeed, the instabil- 
ity of th e bond-centred soliton solutions of the discretised 
Eq.(2.1) is caused by perturbations lying in the conser- 



vative invariant manifold (that is, satisfying the scalar 
nonlinear Schrodinger equation). The instability of the 
bond-centred vector solitons is simply inherited from the 
instability of their scalar counterparts. 



B. Integral considerations 

It is instructive to consider the effect of the perturba- 
tions on the Hamiltonian, momenta an d th e power in- 
tegrals. Substituting Eqs.(|0-(|43]) in ([2^, ([2^ and 



evaluated to be 



{Vu - n)lt^o = Q[{^)H^)dx, (4.10a) 

{Mu-My)\t^Q = -2V2 J q'^{x)(^^dx, (4.10b) 
(Uu - 'l^v)\t^Q = -2V2 J q[{x){a^ + 2(l)^)^dx. (4.10c) 



In Eqs.(4.10) we kept terms only up to the linear order 



in pi^2 and qi^2' 

The eigenvector (4.6) has q[ = q2 = 0; the correspond- 
ing rates (4.10) are all zero. The perturbations associated 



with eigenvectors of this type do not trigger the growth 
or decay of the total power, momentum and the Hamil- 
tonian. They just take the soliton to a nearby solution 
of the scalar Schrodinger equation (3.4); the VT sym- 



metry (the symmetry between the u and v components) 
remains unbroken. 



As for the eigenvectors of the second type, Eqs.( 4.8) 
with nonzero q and p satisfying (4.4b) and (4.9), they 



may set nonvanishing rates of change of the total power, 
momentum and Hamiltonian. Whether the right-hand 



side of Eq.(2.3), (2.5) or (2.6) is zero or not, depends, in 



particular, on the parity of the eigenfunction q{x) of the 



symplectic operator in (4.4b). If (p^q) is an eigenvector 



associated with an unstable eigenvalue (Re /i > 0) and 
suc h th at a particular right-hand side in Eq.(2.3), ( |2.5[ ) 
or (2.6) is nonzero, the corresponding integral in the left- 



hand side will start evolving away from its soliton value. 

In what follows, the perturbations of the integrals per- 
taining to the individual u and v components will also 
prove useful. Restricting ourselves to perturbations as- 
sociated with real eig enval ues / i > 0, and substituting 
Eqs.(|4Tl)-(|0) in and (l2^, we arrive at 



SVu = V2 + 1^ e^' J q[^dx, 
6V. = V2(^^-?j e^'Jq[ 



(pdx, 



and 



Mu = -V^ (^^ + 1) e^* J q'24>.dx, 



My 



1 1 e^W q'2(^xdx 



(4.11a) 
(4.11b) 

(4.12a) 
(4.12b) 



C. Symplectic eigenvalue problem 

Defining X = ax, A = /i/a^, and introducing 

cos 9 



the problem (4.4b) can be written as 



(4.13) 



(4.14) 
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Here Lo,i stand for the scalar Sturm-Liouville operators 



2 sech^X, 
6 sech^X, 



(4.15) 



and we have redenoted qi = g and q2 = f for notational 
convenience. The eigenvalues A and eigenvectors /) 
are generally complex. Positive r] correspond to the soli- 
ton ?/^+ and negative to ?/^_. 

Thus we have reduced a two-parameter stability prob- 
lem to an eigenvalue problem involving a single similar- 
ity parameter. Solitons with different amplitudes and 
in systems with different gain-loss coefficients have the 
same stability properties as long as they share the value 
of 77. (It is fitting to note that a self-similarity of this sort 
was previously encountered in the parametrically driven 
damped nonlinear Schrodinger equation [26^ JTJ. The 
difference of the parametrically driven situation from the 
present setting was that there, the similarity combination 
included two control parameters of the equation whereas 
Eq.(4.13) combines a control parameter with a free am- 



plitude of the soliton.) 

The scalar operators have familiar spectral properties. 
The only discrete eigenvalue of Lq is zero; the associated 
eigenfunction is even: LqZq = 0^ zq = sechX. The opera- 
tor does not have other eigenvalues between and 1 (the 
edge of the continuous spectrum) . The lowest eigenvalue 
of Li is —3; the associated eigenfunction is yo = sech^X. 
The only other eigenvalue is 0; the corresponding eigen- 
function yi = sechXtanhX is odd. 



When 7? = 0, the eigenvalue problem (4.14) coincides 
with the eigenvalue problem for the unperturbed cu- 
bic nonlinear Schrodinger equation. It is important to 
emphasise, however, that the choice 7^ = does not 
correspond to the undamped-undriven situation. (The 
undamped-undriven limit 7 = does not single out any 
particular 77 and is not special in any way.) What the 
value T] = corresponds to, is the turning point 7 = 1. 
As 7 reaches 1 from below, the solitons ?/^+ and tp- merge 
and disappear. 

As in the unperturbed nonlinear Schrodinger, there is a 
four- fold zero eigenvalue at the point r] = 0: A*^^'^'^'^^ = 
0. As 7^ deviates from zero, the four eigenvalues move 
out of the origin. A pertubation calculation (Appendix 
[A| ) shows that two opposite eigenvalues move on to the 
imaginary axis, and the other two move on to the positive 
and negative real axis, respectively: 



1,2) 



±2r?i/2 + 0(7?), A(3.4) = ± A(_^)i/2 + o{v). 

v3 

(4.16) 

The re al an d imaginary parts of eigenvalues of the 
problem (4.14), obtained numerically, are plotted in figj2j 
The two pairs of real and imaginary eigenvalues emerging 
from the origin, are clearly visible. 

Before proceeding to the evolution of the eigenvalues 
as T] grows to large positive respectively negative values, 
it is appropriate to note the position of the continuous 
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FIG. 2. (a) the i magi nary and (b) real part of the eigenvalue A 
of the problem (4.14) as a function of 77. (For the fixed 7, the 



parameter 77 is proportional to the inverse amplitude squared: 
77 — 2 cos 0/a'^.) A real pair of eigenvalues moves onto the 
imaginary axis while the imaginary pair becomes real as 77 
passes through zero (either way). In the 7y > part of the 
figure, the real pair subsequently converges and returns to 
the imaginary axis, so that for 7y > 3, all eigenvalues are pure 
imaginary. In the 77 < part, two complex quadruplets are 
born in two consecutive Hamiltonian Hopf bifurcations. The 
imaginary parts of the quadruplets grow as 77 ^ — cxo while 
the real parts decay. In (a), tinted are the regions Im A > 77 +1 
and Im A < — (77 + 1) filled with the continuous spectrum. 



spectrum the symplectic operator 
-I f Li^T] 



J- 



Lo^T] 



(4.17) 



in (4.14). There are two branches of continuous spec- 
trum, both lying on the imaginary axis of A: A = ±iuj{k)^ 
where uj = 1 -\- r] -\- k'^ . (These are indicated by shading 
in figj2|a).) When 77 > —1, the continuous spectrum has 
a gap7l-^(l + ^),^(1 + ^)]. 
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FIG. 3. The evolution of the perturbed sohton -0+ with the 
amphtude a above the threshold ac- Shown is the magnitude 
of the field u (a) and field v (b). The control parameter is set 
to 7 = 0.1 (for which ac = 0.814); the soliton's amplitude is 
a = 0.820. Note that the spatial interval has been cut down 
to ( — 10, 10) for visual clarity and that only the late stage of 
evolution is shown. In both panels the colour varies from deep 
blue (lowest elevation) to deep red (highest elevation); since 
the maximum values of and |^;| are not equal, the same 
colour corresponds to different values in the two panels. The 
same convention is used in all other double-panel 3D plots. 
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FIG. 4. The evolution of the unstable soliton -0+ may result 
in the formation of a breather. Shown is the magnitude of 
the field u (a) and field v (b). Here 7 = 0.1 and a = 0.9. 



V. HIGH FREQUENCY SOLITON 

In the case of the soliton the eigenvalue problem 
is amenable to simple analysis. The lowest eigenvalue of 
the operator Lq -\- r] equals rj; therefore, when 77 > 0, the 
operator Lq + 77 is p ositive definite and admits an inverse. 
The problem ( |4.14 ) can be written then as a generalised 
eigenvalue problem for a scalar function g{X): 



{Li^r])g = -X^{Lo^vy^g. 



(5.1) 



It is also worth mentioning that the spectrum of 
symplectic operators consists of pairs of opposite pure- 
imaginary values, real pairs and complex quadruplets. If 
A is a real or pure imaginary point of the spectrum, then 
— A is another one; if a complex A is in the spectrum, 
then so are —A, A*, and —A* [28] . 



The operator on the left in (5.1 ) is symmetric, and the 



one on the right is symmetric and positive definite. The 
lowest eigenvalue — A^ in (5.1) is given by the minimum 
of the Rayleigh quotient: 



-A^ = min 



{g\{Lo^v)-'\g)' 



(5.2) 
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FIG. 5. Chart of asymptotic regimes emerging from the 
unstable soUton '0+. The five columns correspond to 7 = 
0.1,0.2,0.5,0.7 and 0.95. The solid line demarcates the 
boundary (5.3); the '0+ soliton is unstable above this line. 



The magenta triangles mark blowups while green circles indi- 
cate the formation of breathers. 



[Here the bra-ket notation is used for the C? scalar prod- 
uct: ^y\z) = J y{X)z{X)dX .] The minimum is positive 
if the lowest eigenvalue of the operator in the numerator 
{v — —3 + T]) is positive: r] > 3. Recalling the definition 



(4.13), we arrive at the stability condition for the soliton 



^4 



a < ar 



r 



(5.3) 



The numerical study of the eigenvalue problem (4.14) 
corroborates these conclusions. As r] grows from to 
positive values, two pairs of opposite eigenvalues, a real 
and a pure imaginary pair, appear from the origin on the 
(ReA,ImA) plane. The real eigenvalues are associated 
with even and imaginary pair with odd eigenfunctions. 
The two pure imaginary eigenvalues diverge to their re- 
spective infinities (figj2|a)). The real pair first grows 
in absolute value, but then the real eigenvalues reverse 
(figj2|b)) and, as r] reaches 3, collide at the origin and 
move on to the imaginary axis. The emerging second 
pair of imaginary eigenvalues diverges to the infinities, 
like the first pair before (fi glg^a )). 

The quantity q{X^t) in ( |472[ ) measures the difference 
between the two components of the vector ([/, V). If the 
eigenvalue problem ( |4.14 ) has a positive eigenvalue A, 
the difference grows monotonically: q = {qi -\- iq2)e^^ ^ 
Therefore, the instability associated with real eigenvalues 
should manifest itself as a monotonically growing asym- 
metry between the two components of the vector field. 
Quantitatively, this asymmetry is measured by the dif- 
feren ce in the increments of the Vu and Vy integrals, 
Eqs.dXTT]). 



Transient and asymptotic solutions emerging as the 
perturbation grows, should inherit the spatial parity of 
the eigenvector (^(X),/(X)) associated with the real 
eigenvalue — that is, should be even in X. 

The imaginary eigenvalues correspond to internal 
modes of the soliton. In Appendix [Bj we derive the 
asymptotes for both pairs of imaginary eigenvalues as 
T] ^ 00: A = ±i{r] + 0.685) + 0{l/r]) and A = 
±i{r]- 1.438) + 0(1/7?). 

Figs [3] and [4] present results of direct numerical simu- 
lations of the ip-^ soliton. For five values of the control 
parameter, 7 = 0.1, 0.2, 0.5, 0.7, and 0.95, we simu- 
lated solitons with amplitudes a above the critical one, 
ttc, given by (5.3). In each of the five cases we have 



identified two possible scenarios of instability growth. In 
one of these, the magnitude of the field u grows without 
bound, while v decreases. (See figjsj) We will be refer- 
ring to this type of evolution as the 'blowup'. Note that 
this asymmetry-growth scenario is in agreement with our 
expectations based on the eigenfunction analysis. 

In the other scenario, the breakup of the unstable soli- 
ton results in the formation of a long-lived oscillatory 
state — a kind of a breather (fig. |4|. Here, the initial 
stage of the evolution is also characterised by the growth 
of asymmetry. Unlike fig|3j it is the v component that 
is growing this time, and u is the one that is decreasing 
(clearly visible in fig. |4|. This "anomalous" growth can- 
not continue indefinitely [this would contradict Eq.(2.3)] 
and eventually the evolution is captured into the breather 
regime. 

Fig|5] sketches the ranges of a which are characterised 
by each of these two types of instability growth. Typi- 
cally, solit ons w ith amplitudes only slightly exceeding the 
threshold (5.3) give rise to breathers whereas the blow- 



up scenario is observed for larger a. The exception from 
this rule occurs for small 7 ~ 0.1 where the blowup and 
breather domains seem to be interlaced in a more com- 
plex fashion. 



VI. LOW FREQUENCY SOLITON 



When ?7 < 0, neither of the scalar operators in (4.14) 
is invert ible on the space. Here, the analysis of the 
eigenvalue problem ( |4.14[ ) has to be done mostly by nu- 
merical means. Our numerical study is summarised in 
the left half of figj2l 

As T] is decreased through zero, a pair of opposite real 
eigenvalues (associated with even eigenfunctions) collides 
and moves on to the imaginary axis. Simultaneously 
(that is, at 7^ = 0) another pair of opposite imaginary 
eigenvalues (also with even eigenfunctions) detaches from 
the continuous spectrum. When 77 reaches the value 771, 



-0.0988, 



(6.1) 



the four imaginary eigenvalues collide, pairwise, and 
leave the imaginary axis forming a complex quadruplet 
A, A*, — A, — A* (with even eigenfunctions). As r] grows to 



larger negative values, the real parts of the quadruplet 
first grow but then start decreasing. At the same time, 
the imaginary parts grow without bound. 

Another pair of eigenvalues colliding as 77 is decreased 
through zero, is pure imaginary. Unlike the colliding real 
pair, the imaginary pair is associated with odd eigen- 
functions. The parity of the eigenfunctions is preserved 
as this pair reappears on the real axis for 77 < 0. As r] 
grows to larger negative values, the absolute values of the 
real A first grow but then start decreasing. 

The existence of the real pair in the interval —l<r]< 
can be established analytically. (See Appendix [C|) It 
might be tempting to expect the pair to converge at the 
origin as r] reaches —1; however the actual bifurcation 
diagram turns out to be more complex (figj2|. In fact 
one can prove that Re A remains nonzero at r] = —1. 
(See Appendix [Cj) 

In the meantime, as r] is approaching —1, the gap in the 
continuous spectrum is shrinking. As 77 reaches —1, the 
two edges of the continuous spectrum meet at the origin 
and the gap closes. Due to the resonance between the 
edge eigenfunctions, a new pair of discrete eigenvalues 
is born at this point. Like the coexisting real pair, this 
pair of opposite real eigenvalues is associated with odd 
eigenfunctions. As r] drops down to 7^2, where 



V2 



-1.168, 



(6.2) 



the newly born pair of real eigenvalues and the real pair 
that has arrived from the origin collide and emerge into 
the complex plane. This is where the second complex 
quadruplet is born. As the negative 77 continues to grow 
in absolute value, the real parts of the eigenvalues making 
up the quadruplet decrease, while the imaginary parts 
grow. 

Thus when 77 < 772, we have two complex quadruplets. 
As T] grows in absolute value, the real parts of the eigen- 
values in both quadruplets decrease whereas imaginary 
parts grow. In the Appendix [B] we derive the asymptotic 
behaviour of the imaginary parts analytically: ImA = 
±{r] - 1.438) + 0(1/77); ImA = ±{r] + 0.685) + 0{l/r]). 
We also show that the d ecrease of the real parts is expo- 
nentially fast. Eq.(4.13) implies then that solitons with 



small amplitude a — even if are unstable — have expo- 
nentially long lifetime. 



r ^ 



exp 



■r 



with some constant <j > 0. 

The nonlinear evolution of the unstable ip- soliton is 
determined by the competition between the eigenvalues 
with positive real parts. According to figj2j there are two 
ranges of the soliton amplitude for each 7: the small and 
the large a. 

For small a [mo re precisely, for a such that |7^| > |7^2| 
with 7^2 as in (6.2)], the soliton has two complex quadru- 
plets in its spectrum — one with an even and the other 
one with an odd eigenvector. Accordingly, the instability 
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FIG. 6. Unstable i/j- solitons with |ry| > 1.168 disperse (a) or 
transform into a long-lived breather (b). In (a), 7 = 0.5 and 
a = 1.075; in (b), 7 = 0.1 and a = 1.152. Both combinations 
of a and 7 correspond to ry = —1.5. Only the ix-component is 
shown in both cases; the evolution of the v field is qualitatively 
similar. 



growth of a small-amplitude i/j- soliton should be ac- 
companied by oscillations. The growth rate of the odd 
perturbation is larger than that of the even one; therefore 
the generic evolution of the instability is expected to be 
dominated by the odd perturbations. 

On the other hand, for very large a — more precisely, 
for 7/1 < 7^ < 0, with 7/1 as in (6.1) — the spectrum has 



only one, real, unstable eigenvalue (with an odd eigen- 
vector). In this range, the instability growth should be 
initially monotonic. The moderately large amplitudes 
(corresponding to 7/2 < t] < rji) are characterised by one 
or two positive real eigenvalues, with odd eigenvector(s), 
and a complex quadruplet whose eigenvector is even in X. 
The odd eigenfunctions have larger growth rates than the 
even ones; hence again, the monotonically growing odd 
perturbations should be dominating the evolution. 
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Thus, depending on the sohton's amphtude, the odd 
perturbations are either the only unstable perturbations 
of the ip- soliton, or its dominant unstable perturbations. 
For these, Eqs.(4.11) give 5Vu = ^Vy = 0. Therefore, the 
odd perturbations do not immediately induce the asym- 
metry of the u and v components of the soliton. On the 
other hand, the momenta Aiu and A4y in Eqs.(4.12) are 
both nonzero, with A4u 7^ Ady. This means that the two 
components, the ''u pulse" and the pulse", will be set 
in motion — they will start moving with unequal veloc- 
ities gradually separating in space. As a result, the "ii 
pulse" will be progressively deprived of the services of its 
power-draining partner, while the pulse" will be cut 
from the power supply by its active counterpart. Even- 
tually, this will set off the growth of the asymmetry in 
the amplitudes of the pulses — the Vu will start growing 
and Vv decre asing. 

Equations (4.12) can tell us whether the emerging 



pulses will move in the opposite or in the same direction 
(yet with different velocities). Indeed, the fragmen ts will 
move in the opposite directions if the momenta (4.12) 
are opposite in signs: 2(7//i) — 1 < 0. Recalling that 
Aa^, this condition can be written as 



2 ^ 27 
a >^. 



(6.3) 



In the complementary region, < 27/A, the u- and v- 
pulses will move in the same direction. 

To verify these predictions, we carried out numerical 
simulations of solitons with a range of amplitudes, for 
small and large values of the gain- loss rate (7 = 0.1, 0.2, 
0.5, 0.7, and 0.95). 

In agreement with the linear analysis, the small- 
amplitude (1 77 1 > 1 7^2 1) unstable solitons were detected 
either to disperse (a process accompanied by the oscilla- 
tion of the soliton profiles, see figj6]^a)) or form long-lived 
oscillatory states (figjGjb)). 

The instability of solitons with large and moderately 
large amplitudes (Ir^l < \r]2\) was seen to grow monoton- 
ically, at least at the initial stage. One of the recorded 
scenarios starts with a spontaneous motion of the two 
components of the soliton in the same direction, with 
slightly different velocities (figj7|), followed by the blowup 
of the ix-component and decay of its v counterpart. This 
behaviour was detected for the amplitudes a lying out- 
side the region ( |6^3[ ). (In particular, the unstable soli- 
ton shown in fig |7] has = 3.12 which is smaller than 
the value 27/A = 4.00 corresponding to 7 = 0.95 and 
T] = -0.2.) 

The other observed evolution starts with the motion 
of the two components in opposite directions. In this 
case the dissociation of the i/j- soliton may be followed 
by a blowup of u and decay of v (figj8| or a formation 
of a pair of breathers (figj9|. Lastly, the linear instabil- 
ity may "just miss" the breathers' basin of attraction in 
which case the nuclei of the two breathers will have their 





FIG. 7. A spontaneous motion of the soliton ip- followed by 
the blowup of its ix-component and decay of v. Here 7 = 0.95 
and a = 1.767 (ry = -0.2). 



tion (6.3). (In particular, figjs^ corresponds to = 9.92 
and 0.42; fig§ to a^= 2.84 and 27/A = 0.31; 

and figliolto = 3.32 and 27/A = 0.31.) 



Our numerical simulations of the unstable regimes are 
summarised in figflTj Triangles mark parameter values 
for which the unstable soliton or the two fragments of 
its break-up were observed to blow up; circles indicate 
simulations that ended in the formation of one or two 
breathers. The solid curve in the figure is the demarca- 
tion line between the small- and large-amplitude ranges 
identified in the line ar a nalysis; the curve is defined by 
T] = r]2 with r]2 as in (6.2). The chart shows a clear corre- 



u component blow up (fig 10). These types of evolution 
were recorded for the amplitudes a satisfying the condi- 



lation between the type of the unstable eigenvalues (real 
vs complex) and the soliton decay product (blowup vs 
breathers). 
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FIG. 8. A dissociation of the soliton '0_ into a pair of pulses 
followed by the blowup of the pulse with large i^- component. 
Shown is the magnitude of the field u (a) and field v (b) . Here 
7 = 0.1 and a = 3.15 (ry = -0.2). 



FIG. 9. A dissociation of the soliton ip- into a pair of 
breathers. Shown is the magnitude of u (a) and contour plot 
of |^;| (b). Here 7 = 0.1 and a = 1.686 (ry = -0.7). 



VII. CONCLUSIONS 

We have examined stability of two families of solitons 
in the two-dimensional VT symmetric coupler with gain 
and loss. The dynamical regimes set off by the instability 
have also been explored. The results of our study can be 
summarised as follows. 

1. Despite the presence of gain and loss, t he b ifur- 
cations occurring in the T^T-symmetric system (2.1) are 
of conservative type. (Linearised eigenvalues cross from 
the imaginary to the real axis, or collide, pair wise, on 
the imaginary axis and emerge into the complex plane.) 
As a result, the soliton instability cannot give rise to lo- 
calised limit cycles (which would be a typical outcome 
of the Hopf bifurcation in a dissipative system). In the 



T^T-symmetric system, the soliton instability either trig- 
gers its blowup (the process where the amplitude grows 
without bound at an exponential rate) or produces finite- 
lifetime breathers. 

2. The soliton stability and internal dynamics are de- 
termined by a single self-similar combination (77) of the 
gain-loss coefficient 7 and the soliton's amplitude a. 

3. The high-frequency solitons with amplitudes smaller 
than ttc, where ac is as in (5.3), are stable, and with 
amplitudes greater than ac, unstable. All low- frequency 
solitons are unstable; however the lifetimes of the solitons 
with small amplitudes are exponentially long so for all 
practical purposes they can be regarded as stable. 

4. The mechanisms of instability of the high- and low- 
frequency soliton are different. The unstable perturba- 
tion of the high-frequency soliton triggers the growth of 



11 





fd 3 



o o 
o o 




0.2 



0.4 



0.6 



0.8 




FIG. 10. A dissociation of the soliton ip- into a pair of pulses 
followed by their blowup. Shown is the magnitude of u (a) and 
modulus of V (b). Here 7 = 0.1 and a = 1.8212 (ry = -0.6). 



asymmetry between the active {u) and lossy (v) com- 
ponents of the soliton, destroying the gain/loss balance 
in the system. The unstable perturbation of the low- 
frequency soliton also upsets the energy balance; however 
this time it is done by splitting the u and v components 
off from their common axis. The difference in the in- 
stability mechanisms is reflected in the difference in the 
products of the high- and low-frequency soliton breakup. 
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Appendix A: Real- imaginary eigenvalue transitions 

The aim of this Appendix is to describe the transitions 
of the pair of eigenvalues of the operator (4.14) from the 



real to imaginary axis and vice versa, as the parameter 
Tj crosses through zero in either direction. 

The numerical analysis indicates that when 77 ap- 
proaches 0, the eigenvalue scales as 77^/^. This suggests 
an expansion of the form 



(Al) 



This expansion will be validated if all coefficient functions 
are found to be bounded and decaying to zero as X ^ 
±00. 



The expansion (Al) is similar to the one appearing in 



the parametrically driven damped nonlin ear Sc hrodinger 
equation [H [27]. The difference of Eq.( |4.14[ ) from the 
eigenvalue problem in [26l [27], is that the parametric 
driving breaks only one of the invariances of the nonlin- 
ear Schrodinger (only the o perator Lq is perturbed). On 
the other hand, in Eq.( |4.14 ), both Lq and Li acquire per- 
turbations. The consequence of this will be the motion 
of two pairs of eigenvalues through the origin on the A 
plane. While one pair is moving from the real to imagi- 
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nary axis, the other pair wih be moving in the opposite 
direction. 



Substituting (Al) in (4.14) and equating coefficients of 



Uke powers of r^, we obtain a chain of equations for A^, 
gn and fn. The order 77^ gives 

^0/0 0, Ligo 0. 

These two equations coincide with equations for the U{1) 
and translational zero modes of the scalar cubic nonhnear 
Schrodinger equation. The bounded solutions are 

/o = Ci sechX, go = C2 sechX tanh X, 

where Ci and C2 are arbitrary constants. 

At the next, 77^^^, order, we get a pair of nonhomoge- 
neous equations 



Lofi — Ai^o, Ligi 
whose solutions are 



-^1/0, 



C C 
fi = ^AiXsechX, gi = ^-^i sechX(l — XtanhX). 

In the context of stability of the scalar cubic nonlinear 
Schrodinger soliton, the generalised eigenvectors (0,/i)^ 
and (^1, 0)-^ would correspond to the Galilean invariance 
of that equation and its soliton frequency variations. 
Finally, the order 77^ gives 



^0/2 = -/o + Ai^i + A2^o, 
Lig2 = —go — ^ifi — ^2/0- 



The solvability condition for Eq. ( A2 ) is 



-(^o|/o)+Ai(zo|^i) =0, 



and the one for (A3) is 



-(^il^o) - Ai(yi|/i 



0. 



(A2) 
(A3) 



(A4) 



(A5) 



In (|A4|)-(|A5|, zq = sechX and yi = sechXtanhX are 
the null eigenvectors of the operator Lq and Li, respec- 
tively. The bra-ket notation is used for the C'^ scalar 
product: 



\z) 



y{X)z{X)dX. 



Substituting for /o,i and ^0,1 in (A4)-(A5), the solv- 
ability conditions reduce to 



A? 



Ci ( ^ - 2 1 = 0, 



Co 



A? 
2 



whence either Ai = ±2 and C2 = 0, or Ai = ±2z/a/3 and 
Ci =0. 



This gives us the the leading-order expressions for the 
two pairs of eigenvalues and eigenvectors of the oper- 
ator (4.14) with small rj. One pair of eigenvalues is 



A = ±2y^ + 0(77); it is associated with even eigenfunc- 
tions: 



/ 



±^/77sechX(l — XtanhX) « x 
sechX ']+0{v). 



This pair moves from the real to imaginary axis as rj 

decreases from positive to negative values. 

The other pair of eigenvalues is A = ±2^^—77/3 + 0(77). 
The corresponding eigenfunctions are odd: 



sechXtanhX ^ ^/ x 



As T] moves from positive to negative, this pair of eigen- 
values translates from the imaginary to real axis. 



Appendix B: Asymptotic eigenvalues as \r]\ -^00 

In this Appendix, we determine the asymptotic be- 
haviour of eigenvalues of the operator (4.14) as 77 tends 
to ±00. 

We expand 



g = go 























V 


^2 



It is convenient to introduce coefBcient functions An{X) 
and B„(X), 

— fn Wn-) ~ {fn H~ Wn) - 



Substituting in (4.14), and equating coefficients of 77 ^, 
gives 



A^i = 1, Ao = 



(Bl) 



at the order 77^, and 

n 

(ii/a - Ao)i5„ = 2sech2x^„ + ^ X^B 

n—m 

(B2) 



ra—1 



2^n+i = 2sech^XSn - (^1/2 + Ao)^n - ^ A^^^- 



for all 71 > 0. Here we have introduced an operator 
Li/2 = -d?/dX'^ + 1 - 4sech^X. 



(B3) 



In Eq.(Bl) we can choose, without loss of generality, 
A_i = 1. [The other root, A_i = —1, defines the opposite 
eigenvalue, —A, with the eigenvector (^, —f)^] 
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Equation (B2) with n = is an eigenvalue problem for 



the operator -L1/2: 

Li/2Bo(X) = \oBo{X). 

The potential — 4sech^X is of the Poschl-Teher variety 
so the eigenvalues can be found exactly. There are two 
discrete eigenvalues, pA and ps' 



PA= OL 



-1.438, pB =3a-4^ 0.685, (B4) 



where a = ^(a/TT — 1). (The corresponding eigenfunc- 
tions are yA = sech^X and ys = sech"~^X tanhX, 
respectively.) Thus the correction Aq may take either of 
these two real values, pA or pB] depending on the choice, 
the function Bo{X) is even or odd. In either case, Bq can 
be chosen real. 

Once A_i, Aq and Bq have been chosen, all higher- 
order coefficients An , Bn and A^ are determined uniquely. 
Furthermore, one can prove by induction that all these 
coefficients are real. The proof proceeds in three steps. 

First, we assume that Bq^ Si, S^, Aq^ Ai^ Ai^ and 
Ao, Ai, A^ with so me i > have been found and are all 
real. Then equation (B3) gives us ^^+1, which does not 



have an imaginary component either. 



Next we turn to the equation (B2 ) with n = £-\-l. The 
solvability condition for Bn is 



^£+l{Bo\Bo) = — ^ Xni{Bo\B£^i-ni) 
m—l 

-2 J Bosech^XAi^idX. (B5) 



This equation expresses A^+i through ^^+1, 
Si, S2, S^, and Ai, A2, A^. By our assumption, 
all these coefficients are real; hence A^+i does not have 
an imaginary part either. Now that the solvability con- 
dition has been satisfied, the nonhomogeneous equation 
(B2) can be solved for S^+i. The right-hand side in (B2) 
includes real .4^+1, So, Si, S^, and Ai, A2, A^+i. 
Therefore, the solution S£+i(X) is real as well. This 
completes the proof. 

Note that if we choose Aq = Pa, all the coefficient 
functions Sn(X) and An{X) are even, whereas if we set 
= pB^ all functions are odd. 

Thus we conclude that in the limit |7^| 00, the real 



part of the eigenvalue in the problem (4.14) is zero to all 
orders in r]~^. This means that the real part is either 



exactly zero or exponentially small in |?7|: ReA ~ e 



-a\r]\ 

a > 0. The imaginary part may assume one of the two 
pairs of values, Im A^ = + Pa + 0{ri~^)] or Im A^ = 
+ Pb + 0(/^~ "^)], where Pa < and p^ > are 
given by Eq.(B4). The eigenfunctions /(X) and g{X) 



associated with the former are even and those pertaining 
to the latter are odd. 



Appendix C: Real eigenvalue in —1 < 77 < 

In this Appendix, we prove the existence of real eigen- 
values of the symplectic operator (4.17) with negative 77, 
and discuss their behaviour as 77 ^ —1. 

If 77 < 0, neither of the scalar operators in (4.14) is 
invertible on the full C'^ space. Fortunately, the sym- 
plectic operator (4.17) is parity preserving and so all its 



eigenfunctions fall into one of the two broad classes: even 
and odd ones. Therefore, when considering eigenvalues A 
with even eigenfunctions, we can restrict the scalar oper- 
ators to the subspace of consisting of even functions. 
When examining the odd eigenfunctions, these operators 
can be restricted to the odd subspace. 

The advantage of the separate treatment of eigenfunc- 
tions with different parity becomes clear when we note 
that the operator + 77 with — 1 < 77 < is positive def- 
inite on the subspace of odd functions — denoted (5 in 
what follows. Therefore, if we confine ourselves to eigen- 
values — A^ of the generalised eigenvalue problem (5.1) 
associated with odd eigenfunctions ^(X), the lowest of 
these "odd" eigenvalue s wi ll be given by the minimum of 
the Rayleigh quotient (|5.2|) on 6: 



-A^ 



mm 



{9\Li+n\9) 



gee {g\{Lo + r])-^\gy 



(CI) 



The eigenfunction yi = sechX tanh X associated with 
the negative eigenvalu e 77 of the operator Li + 77, is odd. 
Hence the quotient in (CI ) can assume negative values in 
6, and its minimum — A^ is negative. We con clude that in 
the parameter region — 1 < 77 < 0, Eq.( 4.14[ ) has a pair of 
opposite real eigenvalues ±A (with odd eigenfunctions). 

Contrary to what one might have expected, this pair of 
eigenvalues does not converge at the origin as 77 reaches 
— 1. The following argument shows that the eigenvalues 
should remain finite as 7^ ^ — 1 + 0. 

First, we note that any function g{X) from (5 can be 
expanded over the complete set of odd eigenfunctions of 
the operator Lq: 



9{X) 



g{k)zk{x)dk. 



(C2) 



Here 



tanhX cos(A:X) + ksm{kX) 



is a continuous spectrum eigenfunction pertaining to the 
eigenvalue E{k) = 1 + 2/c^. The expansion (C2) does not 



include a sum over discrete eigenvalues because Lq has 
only one such eigenvalue and the corresponding eigen- 
function is even. 

Le tting 7/ = — 1 + e, e > 0, and substituting (C2) in 
(CI), the denominator can be written as 



{g\{Lo + v)-'\9) 



\g{k)\^{2e + e)-^dk. (C3) 
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Now consider a subspace & of (3 consisting of func- 
tions g satisfying two conditions. The first require- 
ment is that the function g{X) should decay to zero, 
as X ^ ±00, faster than with any n > 0. This 

condition ensures that the coefficient function Q(k) has 
derivatives of all orders at /c = and hence can be ex- 
panded in a Taylor series centered on that point. 

The second condition is that g{X) should satisfy 



/ 



g{X)zu=o{X)dX = 0, 



(C4) 



where Zk^ 



tanh X is the eigenfunction pertaining 
to the edge of the continuous spectrum of Lq. Owing 
to the constraint (C4), the function Q{k) has a zero at 
= 0; therefore its Taylor expansion has the form Q{k) = 
SiA: + ^2^^ + •••• Consequently, \Q{k)\^k~'^ does not have 
a singularity at /c = 0, the inte gral \ \Q{k)\^k~'^dk is 
convergent and the denominator ( |C3[ ) remains finite as 
e 0. 

If we find a function g in & which renders the numer- 
ator in \C1\ negative, we will obtain a simple bound for 
the eigenvalue A at the point t] = —1\ 



{9\Li-l\g) 
J \g{k)\^k-^dk 



< 0. 



(C5) 



To verify that functions with these properties do exist, 
take, for instance. 



g{X) = sechXtanhX - -XsechX. 



This odd function decay s fa ster than any power of X and 
satisfies the constraint (C4); hence it is in (5 ^ A simple 
integration gives 



(5|(Li-l)|5) = ^-ii^^ 



which is negative (approximately —7.9 x 10 ^). 

This completes our argument. We have shown t hat th e 
positive eigenvalue A of the symplectic operator (4.17), 
which had been proven to exist for —1 < 7^ < 0, does 
not approach zero as 77 approaches —1 from the right. 
(The same is obviously true for the negative counterpart 
of this positive A.) Note that our argument does not rule 
out the existence of pairs of real eigenvalues converging 
at zero as r] approaches —1 from the left. 
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